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Abstract 

We present an ab-initio and analytical study of the Jahn- Teller effect in two diluted magnetic 
semiconductors (DMS) with impurities, namely Mn-doped GaN and Cr-doped ZnS. We show 
that only the combined treatment of Jahn- Teller distortion and strong electron correlation in the 
3d shell may lead to the correct insulating electronic structure. Using the LSDA-|-?7 approach we 
obtain the Jahn- Teller energy gain in reasonable agreement with the available experimental data. 
The ab-initio results are completed by a more phenomenological ligand field theory. 

PACS numbers: 75.50.Pp, 71.23.An, 71.55.-1, 71.70.Ej 
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I. INTRODUCTION 



The Jahn- Teller effect of impurity ions in semiconductors with degenerate ground state 
is well known since long time. Most of its theoretical treatments were based on crystal field 
theory. Actually, the class of diluted magnetic semiconductors regains a lot of interest in 
connection with the search for new materials for spintronics applications. For example, high 
Curie temperatures were predicted in GaN:Mn.- Corresponding experimental^ or ab-initio 
studies^, seemed to confirm these predictions. However, the experimental results are very 
disputed since they might be caused by small inclusions of secondary phases.- And also most 
of the previous ab-initio studies obtained a partially filled band of only one spin direction 
at the Fermi level (half-metallic behavior) which will be shown to be an artifact of those 
calculations. Actually, the Mn ion changes its valence in the chemical series from GaAs:Mn, 
via GaP:Mn to GaN:Mn.~ It is Mn^"*" in GaAs:Mn (for a sufficiently high Mn-concentration)^ 
which leads to hole doping, but it remains Mn^"*" in GaN:Mn. And the Jahn- Teller effect is 
crucial to stabilize Mn^~^. As will be shown below, only the combined treatment of Jahn- 
Teller effect and strong electron correlations leads to the correct electronic structure. The 
electron correlations turn out to be the leading interaction, but the Jahn- Teller effect is 
necessary to break the symmetry. 

We present here a combined ab-initio and crystal field theory of magnetic ions in II- VI or 
III-V semiconductors. As representative examples we treat the Sc?^ ions Mn3+ in GaN and 
Qj.2+ 2nS. The first one is chosen because of its actual interest for spintronic applications 
and the second one since it is a very well studied model system for the Jahn- Teller effect 
of d'^ ions.- Our ab-initio results are in good agreement with the experimental data, but 
only if we include properly the effects of the electron correlations in the 3d shell. For that 
purpose we use the LSDA+f/ method. The resulting electronic structure corrects previous 
electronic structure calculations (which did not take into account the combined effect of 
Jahn- Teller distortion and Coulomb correlation, neither for Mn-doped GaN-'- nor for Cr- 
doped ZuS^iii) in a dramatic way: instead of half-metallic behavior we obtain an insulating 
ground state with a considerable excitation gap. Similar results for GaN:Mn were reported 
earlier but using different methods than in our study.^'^^ In a second step, to make contact 
with the traditional literature on that subject, we connect our ab-initio results with crystal 
field theory. We obtain the complete set of crystal field parameters in good agreement with 
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previous optical measurement s.-ii^*^ 

We treat the host crystals in zinc-blende phase. Both magnetic ions are in the 3d^ 
configuration. The electronic level is split by a cubic crystal field created by the first nearest 
neighbors (ligands) of the transition metal ion. But it remains partially filled and the 
Jahn- Teller effect may occur which induces the splitting of the partially filled level due to 
displacements of the ligands around the transition metal ion. The local symmetry of the 
crystal is reduced and the total energy of the supercell is minimized. The energy gain is 
denoted as Ejt- 



II. METHOD OF CALCULATION 

Our calculations are performed using the full potential local orbital (FPLO)^^ method 
with the LSDA+[/ (local spin density approximation with strong Coulomb interaction)^ 
approximation in the atomic limit scheme. The lattice constants are optimized for the pure 
semiconductors using the LSDA method, = 4.48 A for GaN and ao = 5.32 A for ZnS. 
To study the Jahn- Teller effect, we use a supercell of 64 atoms in zinc blende phase, and 
a 4 X 4 X 4 k-point mesh. The LSDA+?7 parameters are introduced as Slater parameters : 
id), F2 = 7.41 eV, = 4.63 eV and for Cr (see Ref. 



for Mn (see Ref. 



171), = 7.49 eV, 



F^ = 4.68 eV. The F^ parameter is equal to U, and is chosen to be 4 eV. This value is in 
good agreement with other works and it corresponds to the value which gives the maximum 
splitting of the triplet state in D2d symmetry. We have verified that our results are not very 
sensitive to the actual choice of the U parameter. The value U = 4 eV gives representative 
results for a rather large range of U parameters reaching from 3 up to 8 eV. 



FIG. 1: (Color online) Schematic drawing of a tetragonal distortion from (red spheres) 

to (green spheres). 

The transition metal ions (Mn or Cr) substitute one atom in the center of the crystal and 
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there is a complete tetrahedron of Ugands around the magnetic ion. Without distortion, the 
tetrahedron is in the symmetry group (cubic). By symmetry, there are two tetragonal 
(point groups and D2d) and one trigonal {C^^) possible Jahn- Teller distortions. We have 
verified that the most important energy gain is obtained with the pure tetragonal distortion 
where the symmetry reduces from to D2d- That is in agreement with a previous study 
of the Jahn- Teller effect in GaN:Mn where, however, the electron correlation had not been 
taken into account.^ - The schematic displacement of ligands is represented in Fig. [H it is 
defined by 5^ = 5y ^ 5z. 



III. AB-INITIO RESULTS 
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FIG. 2: Partial density of states (DOS) of the two DMS. The results were obtained with 
the LSDA+f/ method (f/ = 4 eV). The dashed line represents the partial "id DOS of the 
transition metal ion and the solid line represents the total DOS. The Fermi level is set to 

eV. 



To study the Jahn- Teller distortion we use a supercell of 64 atoms. We performed a 
series of calculations with different displacements of the ligands around the impurity in 
order to find the configuration which minimizes the total energy of the supercell. The 
preferred configuration is of D2d symmetry and the results for the density of states (DOS) 
are presented in Fig. |2l For Mn-doped GaN, the distance between Mn and N is 1.942 A 
(with cubic symmetry, it is 1.937 A) with 5^ = 1.68 pm and 5^ = 1.76 pm. The Jahn- Teller 
effect induces a lowering of the total energy by 38.13 meV and a splitting of the triplet state 
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by 0.81 eV. For Cr-doped ZnS, the results are similar, the distance Cr-S is 2.314 A (2.304 
A without Jahn- Teller effect), with 5^ = 3.36 pm and 5z = —0.32 pm. The total energy 
decreases by 58.1 meV and the splitting of the triplet level is 1.44 eV. 

Without [/, in the LSDA method, we also find a Jahn- Teller effect for GaN:Mn but the 
energy gain is smaller by nearly two orders of magnitude, only 0.71 meV. In the LSDA 
method all the 3(i levels are located in the gap with a rather small mixing to the 2p orbitals. 
Without Jahn- Teller effect, there is a small band of t2g character which is partially filled and 
one finds half-metallic behavior. The Hubbard correlation opens up a considerable gap. The 
empty state (singlet) is mainly of 3(i character for the two compounds and the doublet, just 
below the Fermi level is mainly of 2p character originating from ligand orbitals. That change 
of orbital character arises since the LSDA+f/ method pushes the occupied "id levels much 
lower in energy than in the LSDA method. As a consequence, the transition to the first 
excited state has a considerable p-d character and should be visible as an optical interband 
transition. That allows a reinterpretation of the optical transition at 1.4 eV (for GaN:Mn) 
which is usually considered as a pure d-d transition.— In agreement with a proposal of Dietl 
it corresponds to a transition from the d'^ configuration to d^ and ligand hole.— 

There are two peaks in the unoccupied, minority DOS at about 2.5 eV for GaN:Mn. 
They correspond to crystal field split ?)d levels which were well seen in X-ray absorption 
spectroscopy. Up to now, they were interpreted by means of a LSDA calculation,— but 
our DOS shows that these peaks occur also in the more realistic LSDA+[/ approach where, 
however, a detailed calculation of the matrix element effects (which was performed in Ref. 



20l ) is still lacking. 

For GaN:Mn, the total magnetic moment is equal to 4/iB, corresponding to S = 2. That 
fits well with the ionicity 3+ for manganese. The local magnetic moment at the manganese 
site is slightly enhanced Mmu = 4.042/1^ which is compensated by small induced magnetic 
moments of opposite sign at the neighboring ligands and further neighbors. 

Another interesting result is presented in Fig. |3l In this case we introduced the local 
Coulomb correlation but no lattice deformation. The local cubic symmetry was however 
broken by the occupation of the triplet state. The calculation still shows an insulating case 
with nearly the same splitting of the triplet state. This observation is not visible within 
the LSDA method. Therefore, we can say that the splitting of the triplet state is mainly 
due to the strong correlation of "id electrons. The Jahn- Teller distortion introduces a small 
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additional splitting of the triplet state : 60 meV for Mn-doped GaN and —11.4 meV for 
Cr-doped ZnS. This interpretation is confirmed by a pure LSDA calculation, because it gives 
an identical value of the splitting with the same ligand coordinates. A negative value of the 
splitting corresponds to a stretching of the tetrahedron along the z axis (Fig. El^b)). 




E [eVl E [eVl 



(a) GaN:Mn (b) ZnS:Cr 

FIG. 3: DOS resulting from the LSDA+f/ calculation (f/ = 4 eV). The dashed line 
represents the DOS without Jahn- Teller distortion, and the solid line represents the DOS 
with Jahn- Teller effect but by breaking the local cubic symmetry. The Fermi level is set to 

eV. 



IV. LIGAND FIELD THEORY 

For a deeper understanding of the Jahn- Teller effect we treat it also in ligand field theory. 
In that theory, the degeneracy of the impurity 3d level is split due to hybridization with the 
neighboring ligands. In the following we neglect the electrostatic contributions which will 
be shown to be justified for GaN:Mn and to a lesser extent for ZnS:Cr. When the 3d ion 
is in the center of an ideal tetrahedron, the cubic crystal field splits this level into a triplet 
state and a doublet state. The distance separating these two levels is denoted as A,. Then, 
the Jahn- Teller effect splits the doublet state into two singlet states, and the triplet state 
into a doublet state and a singlet state. The distance between these levels is denoted as 
AbT^ (shown in Fig. H]). The local Hamiltonian in ligand field or crystal field theory can be 
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expressed as 



CF 



H. 



cub 



tetra 



B^{Ol + hOl) + {Biol + Biol) 



(1) 



where -B^ and are Steven's parameters and Steven's operators, respectively.-^ The first 
part represents the Hamiltonian of an ideal tetrahedron and the second describes the linear 
Jahn- Teller effect. The eigenvalues of Hcf correspond to the Sc? electronic levels of the 
magnetic ion (they represent the spectrum). B^ is the parameter of the cubic crystal field, 
and it is equal to Aq/120. 



'T2 



FIG. 4: Schematic multiplet spectrum for a transition metal ion in df" configuration. From 
left to right, ^D: isolated ion, cubic crystal field splitting, and level splitting due to the 

Jahn- Teller distortion. 



The crystal field Hamiltonian ([T]) has the same form for the one-particle problem (with 



o\ 2132 



parameters i?4, Bl, and Bl) or for the 3(i^ multiplet (with parameters i?4, Bl, and B\ 
The fundamental multiplet state is (L = 2 and S = 2). In this fundamental state, 
the orbital moment for one electron (/ = 2) is equal to the total orbital moment. This 
particular case induces that the one electron spectrum is opposite to the multiplet spectrum 
(ex : i?4 = —B4). In the superposition model22i"— the crystal field Hamiltonian ([T]) can be 
calculated by adding up the hybridization contributions of all the ligands. Then, the matrix 
elements of Hcf with respect to the projection uil of the total orbital momentum L = 2 
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are given by 



% 



(2) 



with 



A , 



:-l)-i5x/4^^224 ^224 



27 



-1)"'^V^^222 r.222 



(3) 



and where C^^^^ are Clebsch- Gordon coefficients and Y^^'^^' are spherical harmonics. The 
i index corresponds to the hgand number. The axis system is defined in Fig. [5l 




FIG. 5: (Color online) Definition of coordinate axes and lattice displacements used in the 

present work. Example of MnN4. 



In identifying the eigenvalues of Hcf with the equivalent values of Kn^,m^, we obtain the 
expressions of the Steven's parameters as a function of the ligand coordinates 



120 



B4 ~ 



84 



12 



12 



(4) 



5° 



-2b2{R)Y,^{e2,e,) 



49 
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with 



= 6 - 30(cos2 62 + cos^ ^3) + 35(cos^ 62 + cos^ ^3) 

= sin"^ 6*2 + sin"^ 6*3 

= -2 + 3 cos^ ^2 + 3 cos^ ^3 



(5) 



9/^4^3 4 

^^^'^^ = 5a;^^^'--3^^-^ 

where 62 and &4(-R) are expressed by means of the Harrison parametrization of "p-d hop- 
ping. The values of r]pda = —2.95 and 77^^,^ = 1-36 are extracted from his book (first 
edition).— In the case of D2d symmetry, the distances i?, between magnetic ion and the 
four hgands are identical. Ap^^ is the charge transfer energy. It is treated as an adjustable 
parameter in this theory. 



TABLE I: Summary of parameters. 



GaN:Mn ZnS:Cr 
our work experiment'^ our work experiment 



As-TT, (meV) 


187 


111 


167.1 


213.9'' iii.e'^ 


EjT (meV) 


59.9 


37 


58.4 


71.3'' , 37.2'= 


Bl (meV) 


-1.88 


-1.05 


-1.26 




Bl (meV) 


-8.24 


-3.98 


-10.15 




A,(eV) 


1.4 


1.37 


0.57 


0.58'' , 0.59'= 


i?(A) 


1.932 




2.303 




(pm) 


1.28 




1.65 




(pm) 


3.79 




3.46 




Aprf (eV) 


2.31 








^2/^4 






3.3 





"optical masurements for wurtzite GaNiMrn Wofos et al, Ref. 



26 



13. 



optical measurements, Vallin et al, Ref. 
'^optical measurements, Kaminska et al, Ref. 
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Hamiltonian ([T]) concerns the electronic degrees of freedom only. The tetragonal distortion 
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results from the coupling to the lattice. As it is shown in Fig. H] any tetragonal distortion 
leads to a splitting of the lowest triplet and to an energy gain. This energy gain is linear 
in the lattice displacements, whereas the vibronic energy loss Eyihronic is a quadratic term. 
One has to minimize the total energy 

2 

\otal = 2 {E{^B2) — E(^E2)) + E^ibronic ■ (6) 

For the sake of simplicity we approximate Embronic by the breathing mode energy, ex- 
tracted from the LSDA+f/ calculation. We find the lattice contribution for GaN:Mn 
to be E^ibronicieV] = 36.7708i?'2 and for ZnS:Cr E^ibronic[ey] = 30.4887i?'2 with R' = 
+ + ^- C)ne should note, that there is no difference between LSDA and 

LSDA+[/ methods for the lattice energy. Also, we have probed the lattice energy for the 
more specific tetragonal mode with no essential difference. The energy gain which is induced 
by the Jahn- Teller effect, corresponds to : 

EjT = min{Atotai) ■ (7) 

The results of the ligand field model are presented in Table [H For GaN:Mn, we have 
exclusively used the ligand hybridization as the microscopic origin for the level splitting. 
The value of Ap^ is adjusted such that Ag equals the experimental value. The results are 
very convincing which proves that the exceptional large value of = 1.4 eV in the case of 
GaN:Mn is dominantly caused by the hybridization energy to the ligands. On the other hand, 
the neglect of the electrostatic contribution is certainly an approximation which shows its 
limits for ZnS:Cr. The procedure described above gives no satisfactory results. We interpret 
this deficiency such that the electrostatic corrections become more important for ZnS:Cr 
which has a much smaller value of Ag indicating a smaller hybridization. As discussed 



in Ref. |27|, the higher order crystal field parameters, and especially are certainly more 
influenced by further reaching neighbors than B4. Therefore, we determine the 64 parameter 
from the experimentally known value of A^ and introduce a second free parameter (62/^4) 
which is fitted to the LSDA+f/ energy gain. The parameter set for ZnS:Cr which was found 
in such a manner (Table H]) is now in good agreement with previous optical measurements. 



10 



TABLE II: Comparison of parameters. When possible, the complete set of parameters was 
calculated from the literatures values by the relations given in Sect V. HSE : 
Heyd-Scuseria-Ernzerhof hybrid functional, GGA : Generalized Gradient Aproximation, 





Gr G : Green-iunct 


ion Galculat 


ion. 








method 


EjT [meV] 


Qe [A] 


V [eV/A] 


huj [cm ^] 


GaN:Mn present work 


ligand field theory 


59.9 


-0.0828 


-1.44 


579 




ab-initio, LSDA+[/ 


38.13 


-0.0562 


-1.35 


680.5 


literature values 


ab-initio, GGA" 


100 


-0.1365 


-1.46 


454 




ab-initio, HSE^ 


184 










experiment, optics'^ 


37 








ZnS:Cr present work 


ligand field theory 


58.4 


-0.0834 


-1.4 


375.8 




ab-initio, LSDA+[/ 


58.4 


-0.0496 


-2.35 


631.5 


literature values 


ab-initio, GFC^ 


185.98 


-0.16 


-2.32 


349.4 




experiment, optics'^ 


37.2 


-0.279 


-0.266 


90 



'^Luo et al. Ref. H 
''Stroppa et al. Ref I 
'^Wolos et al. Ref. L 



iwnrzite GaN:Mn. 



^Oshiyama et al. Ref. ST . 
'^Kaminska et al. Ref. 14. 



V. COMPARISON TO OTHER WORKS 



The lattice displacements obtained from the ligand field theory (Table [T]) fulfill approxi- 
mately the relationship 6^ = 6y = ^6z. That occurs not by accident since this relationship 
characterizes a pure tetragonal mode Qg. In general, up to now, we considered Sx = Sy ^ 6z, 
which corresponds to a mixture of tetragonal and breathing mode Qi, (5^. = 6y = —6z)- 
The normal coordinates of the two modes are defined as Qe = —2-^J^{6x + Sz) and 
Qb = -^{26x — Sz), respectively.— But it is only the tetragonal mode which leads to a 
splitting of the level. That is described in the Hamiltonian which was used by Vallin et 
al ldlj2^ to analyze their data obtained by optical measurements and electron paramagnetic 
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resonance (EPR) 

Hd = VQeee + . (8) 

Here, the V parameter is the Jahn- Teller coupling coefficient, and eg is a diagonal 3*3 matrix 
which describes the splitting of the triplet into the upper ^i?2 doublet and the lower ^i?2 
singlet. Its diagonal elements are 1/2 (corresponding to ^E2) and -1 (corresponding to^i?2)- 
The parameter k describes the lattice stiffness and is connected with the phonon frequency 
by = Kjm where m is the mass of one ligand. Minimizing the energy (|H]) we find the 
relationship V = nQe and the Jahn- Teller energy 

EjT = ^Qe . (9) 

Therefore, if we know the Jahn- Teller energy Ejt and the lattice distortion Qg, we may 
calculate the coupling coefficient V and the phonon frequency u and compare it with other 
works. The comparison is not an ideal one since our LSDA+?7 results do not correspond to 
a pure tetragonal distortion. They contain an important part of the breathing mode which 
comes about due to recharging effects around the impurity being not treated in the ligand 
field theory or in the analysis of the optical data in Refs. Idandy. Nevertheless, we compare 
in Table [III our ab-initio and ligand field results with other theoretical or experimental work 
from the literature. We may remark a rather good agreement for the Jahn- Teller energy Ejt 
between our work and the optical data of Vallin et al for ZnS:Cr.- However, our estimate 
of the phonon frequency is much larger. A detailed analysis of this discrepancy is beyond 
the scope of the present work. We just remark that higher phonon frequencies would mean 
a more profound tendency towards the dynamic Jahn- Teller effect. And there are indeed 
discussions (not for ZnS:Cr, however) whether the experimental data of GaN:Mn should be 
interpreted as a static^^ or a dynamic^ Jahn- Teller effect. Such nonadiabatic effects can, 
however, not at all be treated in density functional based methods. 



VI. DISCUSSION AND CONCLUSION 



Many ab-initio studies of Mn-doped GaN or Cr-doped ZnS, mostly based on the LSDA 
approximation, result in a half-metallic behavior. In those calculations, the Fermi level is 
located in a small 3(i band of majority spin within the middle of the gap of the host semi- 
conductor.-"— We have shown, that such a result is an artifact of the LSDA method and can 
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be repaired by taking into account the strong Coulomb correlation in the 3d shell and the 
Jahn- Teller effect simultaneously. Allowing for a Jahn- Teller distortion only, but neglecting 
the Coulomb correlation leads to a tiny gap at the Fermi level.™ And also the opposite 
procedure of taking into account the Coulomb correlation by the LSDA+f/ method (i.e. 
the same method which we have used) but remaining in cubic symmetry is insufficient.— 
The origin for this deficiency in the given case of a d'^ impurity is a threefold degenerate 
level of t2g symmetry at the Fermi level which is occupied with one electron only. It is 
sufficient to break the local cubic symmetry in the presence of a strong electron correlation 
to obtain insulating behavior with a gap of the order of 1 eV. Such an electronic structure is 
in agreement with the known experimental data for both compounds. Using the LSDA+?7 
method we found the Jahn- Teller energy gain Ejt in good agreement with known optical 
data.i^^ 

In addition to the ab-initio calculations, we developed a ligand field theory to reach a 
deeper understanding of the Jahn- Teller effect. It uses the p-d hybridization between 3d 
impurity and ligands as principal origin of the crystal field splitting.— This hybridization is 
parametrized by the Harrison scheme^ and the lattice energy as obtained from our ab-initio 
calculation is added. The resulting energy gain is very close to the ab-initio results but the 
ligand field theory allows in addition the determination of the complete set of crystal field 
parameters. We find good agreement with the experimental parameter set for GaN:Mn.— 
The comparison is a little bit speculative since the experimental data were obtained for 
Mn in wurtzite GaN. It turns out, that the Mn impurity leads to an additional tetragonal 
Jahn- Teller distortion in addition to the intrinsic trigonal deformation of the host lattice. 
And we find the parameters of this tetragonal distortion close to our results for zinc-blende 
GaN (which can be synthesized, but for which no measurements exist up to now). We have 
observed that our ligand field method works better for GaN:Mn than for ZnS:Cr due to 
the very large cubic crystal field splitting = 1.4 eV in the former case (in contrast to 
Aq = 0.58 eV for ZnS:Cr). 

Our combined ab-initio and analytical study allows also a comparison with the " classical" 
work on ZnS:Cr. The Jahn- Teller efect of that model compound was already studied in the 
seventies in great detail by optical and EPR measurements.-ii^i^^ We obtain an energy gain 
Ejt very close to the experimental data but higher phonon frequencies. 

Finally, we would like to discuss the importance of our results for spintronics applications. 
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First of all, the Jahn-Teller mechanism which we describe is not restricted to the two com- 
pounds of our study. It may occur in all cases where the impurity level is well separated from 
the valence band but only partially filled. The Jahn-Teller effect leads to insulating behavior 
which questions many previous ab-initio studies in the literature. In the case of GaN:Mn all 
available information points to the stability of Mn^+ which hinders an intrinsic hole doping 
by Mn. (Experimentally, one may find Mn^"*" in electron doped samples, which is not well 
suited for spintronics applications either.) One still has the possibility to reach hole doping 
by a second impurity besides Mn (co-doping) . That would allow the classical mechanism for 
ferromagnetism in diluted magnetic semiconductors with S — 2 local moments at the Mn 
sites. 

Our work was supported by a PICS project (No. 4767) and we thank Anatoli Stepanov 
and Audrey Titov for useful discussions. 
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